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COMPARISON OF THERMAL ANALYTIC MODEL WITH EXPERIMENTAL 
TEST RESULTS FOR 30-CENTIMETtR-DIAMETER ENGINEERING 
MODEL MERCURY ION THRUSltR 
by Jon C. Oglebay 
Lewis Research Center 

SUMMARY 

A thermal anal 3 rtic model for a 30- centimeter engineering model mercury-ion 
thruster was developed and calibrated using experimental test results on a pre- 
engineerii^ model 30-centimeter thruster. 

A series of tests were performed at Lewis to simulate a wide range of thermal 
environments on an operating, 30-centimeter engineering model thruster, which was 
instrumented to measure the temperature distribution within it. As a result, both the 
analytic model and the thruster self-heating distribution presented in a previous report 
have been modified. 

This report presents a description of the modified analytic model and comparisons 
of the analytic and experimental results for various operating conditions. Based on the 
comparisons, it is concluded that the anal 5 rtic model can be used as a preliminary 
design tool to predict thruster steady- state temperature distributions for stage and mis- 
sion studies and to define the thermal interfaces between the thruster and other elements 
of a spacecraft. 


INTRODUCTION 

Missions for which the use of large ion thrusters are desired require that the 
thrusters operate over a wide range of power levels and thermal environments (0. 3 to 
5 AU) for as long as 5 years (private communication from C. G. Sauer of the Jet Pro- 
pulsion Laboratory). Tliese missions also require having a large number of thriisters, 
with as few as two or as many as nine thrusters operating at the same time. Because 
of the wide range of power levels and thermal environments, an analytic model is 



required to compute the steady- state temperature distribution within the thruster as 
well as to predict the thermal interaction of the thruster with the environment, adjacent 
thrusters, and the spacecraft itself. 

In addition to the model presented in reference 1, there have been thermal models 
developed for 20-centimeter ion thrusters (refs. 2 and 3) as well as studies on the self- 
heating pattern of operating 15- and 20-centimeter ion thrusters (refs. 4 and 5). In all 
of these studies experimental data were obtained with operating thrusters instrumented 
with thermocouples located in the assumed areas of thermal deposition. The analytic 
models and the thruster self-heating distributions were then adjusted to obtain agree- 
ment of the experimental and predicted temperatures. 

This combined analytic -experimental approach has been used in this study. The 
experimental data were obtained with an operating, 30-centimeter engineering model 
ion thruster (ref. 6). The analytic model presented in reference 1, which was originally 
calibrated using experimental data presented in reference 7, was modified. The self- 
heating distribution was then adjusted to obtain satisfactory agreement with the experi- 
mentally obtained steady-state temperature distributions presented in reference 6. 

Once the self- heating distribution was established, the model was used to compare 
the temperature distributions obtained during simulated multiple thruster operation with 
varied solar flux. Comparisons of the predictions of the analytic model with the data 
obtained during the test program presented in reference 6 will be given. 


THERMAL INVESTIGATION 
Thruster Thermal Characteristics 

The temperature profile within the thruster is a function of the thruster design, the 
operating power level (i.e. , beam current), and the environment in which the thruster 
operates. These factors determine the heat flow paths within the thruster, the amount 
of heat deposited within the thruster, and the sinks to which the heat is rejected. 

Heating of specific thruster components may be required (e.g. , in vaporizers), or, 
may occur due to inefficiencies in the ion beam. The amount of heat supplied to the 
main, cathode, and neutralizer vaporizers (fig. 1) nominally varies from 2 to 7 watts 
each, depending on the beam current (ref. 6). Although this heat causes locally high 
temperatures compared with the main components of the thruster, they are so isolated 
that they have essentially no effect on the temperature distribution within the thruster 
or with the thermal interfaces with the spacecraft or environment. The same applies 
to the cathode and neutralizer tips. 

By far the largest heat source is the discharge power. The amount of heat gener- 
ated in the discharge varies from 111 to 370 watts for beam currents of 1/2 and 2 am- 
peres, respectively (ref. 6). This energy for the most part is deposited in the dis- 
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charge chamber (anode, engine body, etc. ), with the remainder exiting the thruster 
through the apertures in the screen and accelerator grids (fig. 1). This discharge 
power, along with the environment and thruster design, determines the temperature 
distribution within the thruster. 


5 

Thermal Model 

The concept of treating heat flow problems in terms of electrical networks, based 
on the equivalence of the electrical and thermal conductance equations, is well devel- 
oped (ref. 8). The concept basically involves the determination of the thermal resist- 
ance paths within a given network and then numerically determining the resulting tem- 
perature distribution. 

The original thermal network (ref. 1) was reduced from 88 diffusion and 6 boundary 
nodes to 56 diffusion and 6 boundary nodes. A cutaway view of an engineering model 
thruster is shown in figure 1; the nodal layout is shown in figure 2. The nodes and 
materials used are shown in table I. In general, the main components of the thruster 
(i.e., anode, engine body, groimd screen, etc. ) were divided into four circumferential 
quadrants. This approach allows the nodes in each quadrant to see any environment 
desired. This capability was used to simulate multiple thruster operation by surround- 
ing the thruster with variable temperatxire azimuthal shields as will be discussed later. 

The main differences between this model and the original are 

(1) The nodes associated with the neutralizer, cathode vaporizer, cathode, and main 
vaporizer assembles have been omitted. The nodal simulation of these areas was not a 
good representation of the actual hardware; a, more detailed analytic model of these 
areas is required. Fvirthermore, as mentioned before, a previous analysis has shown 
that these assemblies are so isolated from the main components of the thruster that they 
have essentially no impact on the temperature distribution within the thruster. 

(2) The number of nodes on the anode was increased from 8 to 16 (4 circumferential 
by 4 axial) because of the extreme temperature gradient measured during the tests. 

The heat exchange within the thrvister is predominately by radiation. This requires 
that the surface optical properties of the materials be known with a reasonable degree 
of accuracy. However, there can be considerable variation of the optical properties due 
to nonuniform surface finish, the deposition of debris from the operating thruster, and 
any temperature dependence of the surface properties. 

Small samples of the anode and mask were removed, and their optical properties 
measured experimentally between wavelengths of 1 to 15 micrometers. Based on these 
measurements, emissivity values used in the model were 0. 60 for the interior surface 
of the anode, 0. 30 for the external surface of the mask, and 0. 15 for all other thruster 
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surfaces. The vacuum chamber was assumed to have an emissivity of 1. For interior 
nodes radiating through the grids to the cold walls of the chamber, a cavity effect emis- 
sivity of 0.275 was used for the baseplate, and the anode emissivity varied from 0. 675 
to 0.8, depending on the location of the nodes. These values were based on data pre- 
sented in reference 9. 

The radiation exchange factor was approximated«by the following relationship 
developed in reference 8; 



XU XU 

where and are the emissitivities of the i'"" and nodes, A. and A- are the 
1 J th th ^ 

radiating areas of the i and j nodes, and F.. is the geometric configuration factor 

th th 

from the i to the j node . The values of F^j were calculated either by hand or 
numerically using the computer program in reference 10. These inputs were then used 
to calculate the radiation conductor values between the nodes. 

Although radiation is the dominant means of heat transfer within the thruster, the 
linear conduction contribution cannot be overlooked. Therefore, it is necessary that 
the thermal conductivity of the materials as well as the contact conductance across 
joints be known. The thermal conductivities of the materials are well documented in 
the literature and are relatively constant over the temperature range of concern. The 
values for thermal conductivity used herein are presented in table II . Estimating the 
contact conductance across joints is a more difficult task. The conductance between the 
supports connecting the engine body with the anode and ground screen were neglected 
because the contact area is very small compared with the surface area of the nodes. 
Reducing the heat transfer (i. e. , thermal conductivity), to estimate joint conductance 
between other nodes, had essentially no effect on the temperature distribution within 
the thruster. Therefore, joint conductance was neglected in the model. It was also 
assumed that there was no heat transfer between the ground screen and the rear shield 
because of the very poor contact in that area. 

The program chosen to solve the analytical model was the systems improved numer- 
ical differencing analyzer, or SINDA, which is a thermal analyzer program described 
in reference 11. Once the necessary model information was obtained, it was coded in 
SINDA format, and the steady- state temperature distributions of the thruster were 
calculated. Predicting variations in the thermal performance of the thruster due to 
property changes can be done with relative ease. 
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Thruster Self- Heating Distribution 


When the thruster is operating, there are significant heat inputs to various thruster 
components from ion, electron, and plasma radiation heating. The distribution of this 
heat in the thruster is not known precisely and presents another uncertainty in a purely 
anal 3 rtical approach. 

The approach taken here was initially to assume that the heat distribution to the 
various components of the thruster was the same as that used in reference 1 for a 
2-ampere beam current. However, in this series of tests (ref. 6) there was a thermo- 
couple located 1.27 centimeters (1/2 in.) from the grid end of the anode (front) and one 
located 1.27 centimeters (1/2 in. ) from the baseplate end of the anode (rear). The tests 
showed a very sharp temperature gradient from the back to the front of the anode. This 
gradient was undetected in earlier tests because of the lack of instrumentation that close 
to the front of the anode except for one data point (presented in ref. 7) where the mea- 
sured temperature was 655 K. Because of the sharp temperature gradient, the initial 
heating distribution assumed yielded very poor agreement especially in the anode area. 

It was therefore necessary to remodel the anode, as mentioned previously, before 
making any further attempts to estimate the heating distribution within the thruster. 

After the anode was remodelled, the initial heating distribution was adjusted until 
satisfactory agreement between analytic and experimental temperatures was attained for 
a 2-ampere beam cvirrent. The same ratios were then applied for beam currents of 1 
and 1/2 ampere. The final values and percentages of the total heat input to the various 
components for the three beam currents are shown in table HI along with the total dis- 
charge power. The total heat input is less than the total discharge power because some 
heat escapes through the apertures in the grids. 

As can be seen, the majority of the heat goes to the front of the anode (59 percent); 
whereas in references 1 and 2, only 18 percent was assumed to reach the front of the 
anode. The values for the other components of the thruster are as much as a factor of 
5 less than the distributions presented in references 1 and 2. Based on the above com- 
parisons, it appears that the thrusters are not heated as uniformly as has been pre- 
sented in the past. It is felt, however, that the approximations in table H can be used 
for other size thrusters of the divergent magnetic field type. 


RESULTS AND DISCUSSION 
Calibration Tests 


The anal 5 dic model was calibrated using data for a 2 -ampere beam current for a 
single thruster in a vacuum chamber surrounded by Liquid- nitrogen- cold walls both with 
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and without simulated solar flux. To be acceptable, the analytic model temperatures 
were to agree to within ±10 percent of the experimental results. This level of accuracy 
was chosen because the temperatures are not near the material limits and, therefore, 
better accuracy is not required for the safety of the thruster. Also, this level of ac- 
curacy characterizes the thruster as an element of a spacecraft. Once the heating dis- 
tribution shown in table in was established for 2 amperes and no solar incidence, runs 
were made using the distributions for 1- and 1/2-ampere beam currents. These beam 
currents represent full, half, and quarter power. 

The comparison of experimental and analytic temperature distributions for these 
three beam cvirrents are shown in table IV. Shown on this and ensuing tables are the 
experimental and analytic temperatures, the percent error between the two, and the 
average of the absolute values of the errors not including the rear shield. All of the 
analytic predictions are well within 10 percent of the experimental results, with the 
average absolute error ranging from 2. 64 percent at 2 amperes to 4. 26 percent at 
1/2 ampere. Generally, the calculated values are less than the measured values. The 
largest errors are generally in the grovind screen temperatures at a 2-ampere beam 
current. As can be seen, the 2-ampere-beam- current model overpredicts the tempera- 
ture by as much as 29 K, the 1-ampere beam current model is within the acciiracy 
limits, and the 1/2-ampere beam current model underpredicts the temperature by as 
much as 14 K. These differences could be due to uncertainties in estimating view fac- 
tors and effective radiating areas for curved, perforated surfaces, difficulties in meas- 
uring accurate temperatures, or some other phenomena. In any case, the grotmd 
screen temperatures are at least 100 K lower than the other components, and the errors 
are tolerable. As will be shown later, when the thruster is enclosed by a shield and 
the temperatures are much warmer, the differences sharply decrease. 

The error in the rear of the engine body and the anode temperatures increase with 
decreasing beam current, indicating that perhaps the self-heating distribution is not 
directly proportional to the discharge power. The model consistently underpredicts the 
mask temperature. The reason for this is probably that the emissivity could vary sub- 
stantially aroimd the mask, as was indicated by visual inspection of the thruster. (The 
emissivity of the mask sample measured in one location was approximately 0. 3. ) 

The rear shield temperatures were set to the average experimental values because 
of the complex nature of the support structure and instrumentation leads, etc. (refs. 6 
and 7). Attempts were made to estimate the heat rejection from the rear shield. How- 
ever, the model experienced convergence problems for most cases, and the effort was 
abandoned. For the cases that converged an apparent rear shield emissivity of approx- 
imately 0.40 yielded satisfactory results. However, in a more realistic situation, the 
environment the rear shield would view probably would not be as cluttered with support 
structvire, instrumentation leads, etc. , and the rear shield interaction with that en- 
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vironment could be incorporated in the model. 

. In order to establish the effective solar absorptivity of those components of the 
thruster that were subjected to solar incidence (i. e., mask, accelerator grid, and 
baseplate), the model was calibrated with a 2-ampere beam current and the same self- 
heating distribution as before but with a solar incidence of 2 suns axially on the grid end 
(front) of the thruster. Initially, the values of solar, absorptivity used were typical 
handbook values for the materials. These values were later adjusted to obtain satis- 
factory results. It was assumed that half of the solar flux incident on the grid area went 
through the holes to the baseplate because the grids open area was 50 percent. Table V 
shows the solar heating distribution to the thruster and the assumed values of solar 
absorptivity used in these studies. The values of the solar absorptivities are all within 
25 percent of the initial handbook quoted values used for the thruster materials. 

A comparison of the experimental and analytic temperatures for the three beam 
currents are shown in table VI. As in the cases with no solar flux, the analytic predic- 
tions are well within 10 percent of the experimental results, except for the ground 
screen temperature, although they were a little higher than the no solar flux cases. The 
average absolute error ranges from 4. 16 percent at 2 amperes to 6. 22 percent at 
1/2 ampere. The model consistently overpredicts the ground screen temperatures. 

The probable reason for this (in addition to the other possibilities mentioned earlier) is 
that the assumed conductive heat transfer between the ground screen and the mask is 
probably too high. No attempt was made to determine the sensitivity of the ground 
screen temperature to the conductive heat-transfer value because the disagreement, 
again, sharply decreases when a shield is aroimd the thruster. 

Comparing tables IV and VI, the largest increase in temperature due to the solar 
flux occurs on the mask, which is a solid area directly facing the solar source. The 
analytical increase in mask temperature is 20 to 30 K greater than the experimental in- 
crease, depending on the beam current. As the beam current decreased, the increase 
in temperature became greater not only for the mask, but for the other components of 
the thruster as well. Omitting the ground screen and mask temperatures, the average 
experimental increases in component temperature were 30, 34, and 53 K for beam cur- 
rents of 2, 1, and 1/2 ampere, respectively, while the analytic increases were 23, 29, 
and 49 K. The solar flux was more in evidence at the lower beam currents because the 
thruster was cooler due to the lower discharge power. 


Multiple Thruster Simulation 

After the model was calibrated, it was compared with an experimental simulation of 
multiple thruster operation. Multiple thruster arrays were simulated by surrounding the 


7 



thruster with two different variable temperature azimuthal shields. The first configura- 
tion had the thruster completely surrounded by a shield as shown in figure 3; the second 
configuration had the thruster surrounded by a 270° shield with the neutralizer qviadrant 
open to the vacuum chamber walls (fig. 4). These shields allowed the thruster to radi- 
ate primarily to a warm body rather than the liquid-nitrogen-cold walls. This would 
permit a thermal evaluation of an operating thruster surrounded by either operating or 
nonoperating thrusters. The shields were instrumented with thermocouples to monitor 
the shield temperatures. Strip heaters were attached imiformly to the bottom of the 
shields to control the temperature. Samples of the heater were removed, and optical 
property measurements made on the surface facii^ the thruster. The resultant emis- 
sivity was approximately 0. 60. This value was used in the analysis. 

The thruster was operated at the same beam cvirrents as in the calibration tests, 
both with no solar incidence and with up to 2 suns of solar incidence. A series of 
15 tests, 5 each at the three beam currents, were compared with analytic predictions. 

A comparison of the results for the 2-ampere beam current are shown in table VII. 

All of the analytic temperatures are well within 10 percent of the experimental values, 
including the ground screen, which indicates that the interaction between the ground 
screen and chamber environment is not being treated correctly without the shields 
around the thruster. In fact, the errors in the analytic predictions for the most part 
are lower than for the calibration cases. The maximum average absolute error is 
1. 55 percent, with the largest error in these cases occurring on the front of the anode. 

The only major differences between these shielded cases and the calibration cases 
are the much better ground screen agreement and the higher (due to the presence of the 
shields) thruster temperatures. The temperature around the thruster, excluding the 
grid end, was increased from liquid -nitrogen temperature (80 K) to nearly the thruster 
temperatvu-e (500 K). This essentially eliminated the heat transfer from the sides of 
the thruster to the environment. This 420 K increase caused a maximum average ex- 
perimental thruster temperature rise, excluding the ground screen, the front of the 
anode, and the rear shield, of only 68 K in the 360° shield case with 2 suns. The corre- 
sponding anal 3 dic temperature rise was 59 K. This means that the main heat-rejection 
mechanism is by radiation out of the accelerator end of the thruster to the cold walls of 
the chamber . 

Similar results for the 1- and 1/2-ampere beam currents are shown in tables vni 
and K, respectively. As in the calibration tests, the errors generally increase with 
decreasing beam current but still remain well within 10 percent. 

The temperature comparisons of the experimental and analytic results for the 
2-ampere beam current are also given in figures 5 to 11. In these figures the average 
temperatures of the engine body and ground screen presented in tables IV and VI and VII 
are shown. 
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A computer listing of the SINDA input, execution procedure, and output for a typical 
run with a 2-ampere beam cimrent is presented in the appendix. It should be noted that 
all input values are in English imits because SINDA requires that the input temperatures 
be in degrees Fahrenheit for problems involving radiation. 

CONCLUSIONS 

An analytic thermal network for a 30-centimeter engineering model mercury-ion 
thruster has been modified and calibrated against experimental data obtained on an op- 
erating 30-centimeter engineering model thruster. Experimental data were generated 
at full, half, and quarter beam currents over a wide range of thermal environments. 
Comparisons of the anal 5 rtic and experimental temperatures show that the analytic model 
agrees well, within 10 percent of the experimental results (about 4 percent on the aver- 
age), for the main components of the thruster. 

From the results of this study, it was concluded that 

1. A realistic thermal model has been constructed to represent the thruster tem- 
perature distribution with an expected accuracy of 10 percent. 

2. Heat losses to^ various components of the thruster from the plasma are much 

I ' 

different from those previously estimated. 

3. The model can be used to analyze interactions between thruster arrays, thermal 
environment, and the spacecraft itself. 

4. The model can be used for preliminary design, flight performance, and other 
applications . 

5. The main heat rejection mechanism is by radiation out of the accelerator end of 
the thruster. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 16, 1977, 

506-22. 
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APPENDIX - SAMPLE OF SINDA INPUT, EXECUTION PROCEDURE, AND OUTPUT 


SINDA Input 


BCD 3THERMAL SPCS 

BCD 9 30 CM. EMT MODEL (NO HEATER) 

END 

BCD 3N0DE DATA 
PEM 

REM NODE NUMBER, INITIAL TEMPERATURE (F), AND NODE 

REM THERMAL CAPACITY (BTU/F), NEGLECTED IN STEADY STATE 

REM 

101. 83. . .0018, 102 ,80. ,..0016, 10 3 , 80 . , . 00 1 8 , 1 04 , 8 0 . , . 00 1 8 
105 ,8D. , .0023, 106,80. ,.0023, 107, 80 ., .0023 , 1 08 , 80 . , . 002 3 

109.80. . .0023, 1 ID ,80. , .0023, 1 1 1 , 80 . , . 002 3 , 1 12, 80., .0023 

113,80 . ,.0018, 1 14 ,80. , .0018, 1 1 5 , 80 . , . DO 1 8 , 1 1 6 , 8 0 . , . 00 1 8 
9,eC.,.G609, 10, 80. , .0609, 1 1 , 80 . , . 060 9 , 12,80. ,.0609 

13.80. . .0100, 14 ,80. , . 0100 , 1 5 , 80 . , . 01 00 , 16,80. ,.0100 

17 .80 . . .01 1 6, IB, 80. ,.0116, 1 9 , 8 0 . , . 0 1 1 6 , 20,80. ,.0116 

21 .80. . .03024, 22, 80. , .03024 , 23 , 80 . , . 0 3024 , 2 4 , 80 . , . 0 3024 

25.80. . .05316, 2 6 , 80 . , . 0 53 16 , 27 , 30 . , . 05 3 16 , 28 ,80 ., .05316 

29.80. . .0255, 30 , 8 0 . , . 02 55 , 31 , 8 0 . , . 02 55 , 32, 80. ,.0255 

33 .80. , .01 32, 34 , 80. ,.01 32 , 35 ,80. , .01 32 , 36,80. ,.0132 

37.80. . .0048, 38 ,80. , .0048 , 39 , 8 0 . , . 00 48 , 40,80. ,.0048 
PEM REAR SHIELD BOUNDARY NODE (F) 

"■4 6,211.,0. 

74 .80. , .042 5, 75 ,80. , .0425 , 76 , 80 . , . 04 25 , 77,80. ,.0425 

85 . 80. . .0107, 86 ,80. , .0107 , 87 , 8 0 . , . 0 1 07 , 88,80. ,.0107 
REM ENVIRONMENTAL BOUNDARY NODES (F) 

-501 ,-310. ,0. , -502,-310. ,0., -503 , - 3 1 0 . , 0 . , - 5 04 , - 31 0 . , 0 . 
-505,-310. ,0. 

END 


OiLATIVt ><ODf VUHSeSS ACTUAL N00£ NUnBfRS 


1 

THSU 

13 

101 

102 

103 • 

104 

1C5 

1D6 

107 

lOB 

109 

110 

1 1 

THRU 

23 

1 1 1 

112 

113 

114 

ns 

116 

9 

10 

11 

12 

21 

THRU 

33 

1 3 

14 

15 

16 

17 

18 

1 9 

20 

21 

22 

71 

THRU 

43 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

4 1 

THRU 

63 

33 

34 

35 

36 

37 

38 

39 

40 

74 

75 

5 1 

THRU 

63 

76 

77 

85 

86 

87 

88 

46 

501 

502 

5C3 

bl 

THRU 

62 

ST*) 

5C6 









MOOT 

analysis. . 

. . DIFFUSION r 

56, 

ARITHMETIC t 

0, 

BOUNOAFY r 

6 • 

TOTAL r 

62 




BCD 3CONDUCTOR DATA 
REM 

REM LINEAR CONDUCTOR NUMBER, CONNECTING NODES AND VALUE (BTU/HR.F) 
REM 

30 0.1 01, 105, 10 2, 106, 103, 107, 104, 108, 109, 113, 110, 114, 111, 115 

112.116. . 0322 

30 1.105. 109. 106. 11 0.107. 111. 108. 112. . 0161 

3001,101,102,102,103,103,104,104,101,113,114,114,115,115,116 

11 6 . 1 1 3. . 0004 

3002,105,106,106,107,107,108,108,105,109,110,110,111,111,112 
11 2, ID 9,. 0 0 06 

302,9, 10,1 0,11 ,11,12,12,9, .0 092 8 

303. 1 3, 14, 14, 15, 15, 16, 16, 13, .00742 

30 5. 9. 13. 10. 14. 11 . 15. 12. 16 . . 0399 
30 6, 9, 74, ID, 75, 11, 76, 12, 77, .0377 
30 7, 13, 17, 14, 18, 15, 19, 16, 20, .1795 
30 8, 17 , 21, 1 8, 22, 19, 23, 20, 24,. 141 
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30 9, 17, 18, 18, 19, 19, 20, 2D, 17, .01095 

31 0.21. 22. 22. 23. 23. 29. 24. 21.. 068 

311. 25. 26. 26. 27. 27. 28. 28. 25.. 0 74 

31 2. 29. 30. 30. 31. 31. 32. 32. 29.. 0184 

313. 33. 34. 34. 35. 35. 36. 36. 33.. 0372 

314. 33. 37. 34. 38. 35. 39. 36. 40. . 29 

31 5. 37. 30. 38. 39. 39. 40. 40. 37.. 009 

336. 74. 76. 75. 76. 76. 77. 77. 74. . 0 04 6 

REM GROUND SCREEN DISCONNECTED FRM RE4R SHIELD 

347.33.46.34.46.35.46.36.46.0. 

34 8, 9, 29, ID, 30, 11, 31, 12, 32, .0145 

349. 37. 86. 3 8.36. 39. 87.40.88. . 106 

350. 85. 66. 86. 87. 87. 88. 8 8. 65. . 0053 

REM 

REM INTERNAL RADIATION CONDUCTOR NUMBER, CONNECTING NODES, 

REM VALUE INDICATOR, 1.- CONNECTED, Q.- NOT CONNECTED 

REM ACTUAL VALUE CALCULATED LATER 

REM 

-401, 101,102,102, 103,103, 104, 104,101, 113, 114, 114, 115, 115, 116 

116.113.1. 

-402,105,106,106,107,107,108,108,105,109,110,110,111,111,112 

112.109.1. 

-4 03, 101, 103, 102, 104, 11 3, 115, 114, 116,1. 
-404,105,107,106,108,109,111,110,112,1. 

-405,101,105,109,113,102,106,110,114,103,107,111,115,104,106 

112.116.1. 

-406,105,109,106,110,107,111,108,112,1. 

-407,101,106,102,107,103,108,104,105,113,110,114,111,115,112 

116,109,102,135,103,106,104,107,101,108,114,109,115,110 

U6,lll,113ill2,l. 

-4 08, 101, 107, 102, 108, 103, 105, 1 04, 106, 113, 111, 114, 112, 115fl09 

116,110,1. 

-4 09, 105, no, 106, 111, 107, 112, 108, 109, 106, 109, 107, no, 108, 111 

105.112.1. 

-410,105,111,106,112,107,109,108,110,1. 

-4 11, 101, 109, 1 02, no, 103, 111, 104, 112, 113, 105, 114, 106, 115, 107 

116.108.1. 

-4 12,101,1 1 0,1 02,1 1 1,103,1 12,104-, 109,1 1 3,106,1 14,107,11 5,108 
116, 1 0 5, 101, 112, 102, 109, 103, 110, 104, 111, 113, 108, 114, 105 

115.106.116.107.1. 

-413,101,111,102,112,103,109,104,110,113,107,114,108,115,105 

116 . 106 . 1 . 

-4 14, 101, 113, 102, 114, 103, 115, 104, 116,1. 

-415,101,114,102,115,103,116,104,113,101,116,102,113,103,114 

104,115,1. 

-4 16, 101, 115, 102, 116, 10 3, 113, 104, 114,1. 

-4 17, 101,17, 102, 18, 103, 19, 104, 20,1. 

-4 16, 105, 17, 10 6, 18, 107, 19, 108,20,1. 

-4 IV, 109, 17, no, 18, 111, 19, 112, 20,1. 
-420,113,17,114,18,115,19,116,20,1. 

-4 21, ID 1,1 8 , 131, 20, 102, 17, 102, 19, 103, 16, 103, 20, 104, 17, 104, 19,1. 
-4 22, 105, 1 8 , 105, 20, 106, 17, 106, 19, 107, 18, 107, 20, 108, 17, 108, 19,1. 
-4 23, 1 09, IE, 109, 20, no, 17, no, 19, 111, 18, 111, 2 0,112, 17, 112, 19,1. 
-4 24, 113, 16, 113, 20, 114, 17, 114, 19, 1 15, 18, 11 5, 23, 116, 17, 116, 19,1. 
-425,101,19,102,20,103,17,104,18,1. 

-4 2b, 105, 19, 106, 20, 107, 17, 108, 18,1 . 

-4 27, 109, 15, 110, 20, 111, 17, 112, 18,1. 

-4 26, 1 13, n, 114, 20, 115, 17, 116,18,1. 

-4 29, 101,21 ,102, 22, 103, 23, 104, 24,1 . 
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-4 3D, 1 05, ?.l ,136,22,107,23, . 

-4 31, 109,21, 110, 22, 111, 23, 112, 24,1. 

-4 32, 113, 21, 114, 22, 115, 23, 116, 24,1. 

-4 33, 101, 22, 101, 24, 102, 21, 102, 23, 103, 22, 103,24, 104, 21, 104, 23,1 
-4 34, 105, 22, 105, 24, 106, 21, 10 6, 23, 107, 22, 107, 24, lOy, 21, 108, 23,1 
-435,109,22,109,24,110,21,110,23,111,22,111,24,112,21,112,23,1 
-4 36, 113, 22, 113, 24, 114, 21, 114, 23, 115, 22, 115, 24, 116, 21, 116, 23,1 
-437,101,23,102,24,103,21,104,22,1. 
-438,105,23,106,24,107,21,108,22,1. 

-4 39, 1 09, 23, no, 24, 111, 21, 112, 22,1. 
-440,113,23,114,24,115,21,116,22,1. 

-4 41, 101, 29, 102, 30, 103, 31, 104, 32,1. 

-4 42, 105, 29 , 106, 30, 107, 31, 108, 32,1. 

-443,109,29,110,30,111,31, 112,32,1. 

-4 44, 113, 29, 114, 30, 115, 31, 116, 32,1. 

-445, 101, 30, 101, 32, 102, 29, 102, 31, 103, 30, 103, 32 , 104,31,104,29, 
-4 4b, 105, 30, 105, 32, 106, 29, 106, 31, 107, 30, 107, 32, 10c, 31, 108, 29,1 
-447,109,30,109,32,110,29,110,31,111,30,111,32,112,31,112,29,1 
-4 48, 113, 30, 113, 32, 114, 29, 114, 31, 115, 30, 11 5, 32, 116, 31, 116, 29,1 
-4 49, 101, 31 , 102, 32 , 103, 29, 104, 30,1. 

-4 5 0,105, 31 , 106, 32, 107, 29, 108,30,1 . 

-4 51, 109, 31 , 110, 32 , 111, 29, 112, 30,1. 

-4 5 2, 113, 31, 114, 32, 115, 29, 116, 30,1. 

-461, 101, 9, 102, 10, 103, 11, 104, 12,1. 

-4 62,1 05,9, 105, ID, 107, 11, 108, 12,1. 
-463,109,9,110,10,111,11,112,12,1. 

-4 64, 113, 9, 114, 10, 115, 11, 116, 12,1. 

-4 65, 101, 10, 101, 12, 102, 9, 102, 11, 103, 10, 103, 12, 104, 9, 104, 1 1 ,1. 
-4 66, 1 06, 1 0 , 105, 12 , 106, 9, 106, 11 , 107, 10, 107, 12, 108, 9, 108, 1 1 , 1 . 
-4 67, 109, 10 , 109, 12, no, 9, no, 11, 111, 10, 111, 12, 112, 9, 112, 11,1. 
-4 68, 113, 10 , 11 3, 12 , 114, 9, 1 14, 11, 115, 10, 116, 12, 116, 9, 116, 11 ,1. 
-469,113,13,114,14,115,15,116,16,1. 

-4 70, 109, 1 3 , no, 14, in, 15, 112, 16,1. 

-471,113,14,113,16,114,13,114,15,115,14,115,16,116,13,116,15,1 
-4 72, 109, 14, 109, 16, no, 13, no, 15, 111, 14, 111, 16, 112, 13, 112, 15,1 
-1029, 17, 29, 13, 30, 19, 31, 20, 32,1. 

-10 30, 17, 31, 18, 32, 19, 29, 20, 30,1. 

-1031, 17, 30, 17, 32, 18, 29, 18, 31, 19, 3D, 19, 32, 20, 31, 20, 29,1. 

-1034, 17, 21 , 18, 22, 19, 23, 20, 24,1 . 

-ID 35, 17, 22, 17, 24, 18, 21, 18, 2 3, 19, 22, 19, 24, 20, 23, 20, 21,1. 

-10 36, 17, 23, 18 , 24, 19, 2 1,20, 22,1 . 

-1037, 17,18,18,19,19,20,20,17,1. 

-1038,17,19,16,20,1. 

-ID 39, 21, 29 , 22, 30, 23, 31, 24, 32,1 . 

-ID 40, 21, 30, 21, 32, 22, 29, 22, 31, 23, 30, 23, 32, 24, 31, 24, 29,1. 
-1041,21,31,22,32,23,29,24,30,1. 

-1045, 1 3, 85, 14, 86, 15, 87, 16, 88,1 . 

-1046, 25, 86 , 25 , 36, 26, 85, 26, 67, 27, 86, 27, 88, 28, 87, 28, 65,0 . 

-1047, 21, 25 , 22, 26, 23, 27, 24 , 29,1. 
-1052,17,13,18,14,19,15,20,16,1. 

-1053,17,14,17,16,18,13,18,15,19,14,19,16,20,13,20,15,1. 

-10 54 , 9,33 , 10, 34 ,1 1 ,35 , 12, 36,1 . 

-1055, 13, 37, 14, 38, 15, 39, 16, 40,1. 

-1056, 13, 9, 14, 10, 15, 11, 16, 12,1. 

-1063,29,46,30,46,31,46,32,46,1. 

-1122,74,13,75,14,76,15,77,16,1. 

-1 123, 74, 14, 74, 16, 75, 13, 75, 15, 76, 14, 76, 16, 77, 15, 77, 13,1. 

-1124,74,37,75,36,76,39,77,40,1. 

-1 125, 74, 38 , 74, 40, 75, 37, 75, 39, 76, 36, 76, 40, 77, 39, 77, 37,1 . 
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-1143,13,40,15,40,1. 

-1144,9,36,11,35,1. 

-1 145, 13, 38, 14 , 37, 14, 39, 15, 38, 16, 37, 16, 39,1. 

-1146, 9, 34, 10, 33, 10, 35, 11, 34, 12, 33, 12, 35,1. 

REM 

REM EXTERNAL RADIATION CONDUCTOR NUMBER, CONNECTING NODES, 

REM VALUE INDICATOR 

REM 

-20UD, 10 1,5 05, ID 2, 505, 10 3, 505, 10 4, 50 5,1 , 

-2001, 105, 5 05, 106, 505, 107, 505, 108, 505,1. 

-20011,109, 505,110,505,111,505,112,505,1 . 
-20012,113,505,114,505,115,505,116,505,1. 

-2002,9,501 ,1. 

-2003,10,502,1 . 

-2004,11,503,1 . 

, -2 0 b:5 ,12,504,1 . 

-2006,13,501,1 . 

-2007,14,502,1 . 

-2008, 15,503,1 . 

-2009,16,504,1 . 

-2010,17,505,18,505,19,505,20,505,1. 

-2011,21,505,22,505,23,505,24,505,0. 

-2012,25,505,26,505,27,505,28,505,1. 

-2013,29,505,30,505,31,505,32,505,1. 

-2014,33,501,1 . 

-2015 , 34,502,1. 

-2016,35,503,1 . 

-2017,36,504,1 . 

-2018,37,501,1 . 

-2019,38,502,1 . 

-2020, 39,503,1 . 

-2021,40,504,1 . 

-2041,85,505,86,505,87,505,88,505,1. 

-2043, 33,502, 1 . , -2044,33,504,1., -2045,34,501,1. 
-2046,34,503,1., -204 7 , 3 5 , 50 ? , 1 . , -2048,35,504,1. 
-2049,36,503,1., -2050,36,501,1., -2051,37,502,1. 
-2052,37,504,1., -2053,38,501,1., -2054,38,503,1. 
-2055,39,502,1., -2056,39,504,1., -2057,40,503,1. 
-2059,40,501,1 . 

-2059,33,505,34,505,35,505,36,505,1. 

-2 065, 9, 50 5 , 13, 505, 11 , '^05, 12,505,1 . 
-2066,13,505,14,505,15,505,16,505,1. 

-2 067, 74, 505, 75, 50 5, 76, 505, 77, 505,1. 
-2070,37,505,36,505,39,505,40,505,1. 
-2071,85,501,86,502,87,503,0. 

-2072,88,504,3. 

-2073,74,501,75,502,76,503,1. 

-2074 ,77,504,1. 

END 

OELATIVE conductor "JUHBEPS actual COMOUCTOP N'UrPtPS 


1 

THRU 

1 D 

3ra 

301 

3001 

3002 

3C2 

303 

505 

306 

3:^7 

zna 

1 1 

THRU 

23 

3C9 

310 

311 

312 

313 

314 

315 

336 

347 

343 

21 

THRU 

33 

349 

350 

401 

402 

403 

404 

4C5 

406 

• 407 

405 

31 

THRU 

43 

409 

413 

411 

4 12 

415 

414 

415 

4 1 b 

4 1 7 

4 1 5 

tt 1 

THRU 

53 

419 

420 

421 

422 

423 

«24 

425 

426 

427 

4?r. 

5 1 

THRU 

63' 

429 

4 30 

4 31 

432 

433 

4 34 

4 35 

436 

437 

4 33 

61 

THRU 

73 

4 39 

440 

441 

442 

44 3 

444 

445 

446 

447 

445 

71 

THRU 

83 

449 

450 

451 

452 

461 

462 

465 

4 64 

46b 

46t 

S 1 

THRU 

93 

467 

46b 

469 

4 70 

471 

4 72 

1029 

1030 

IL’l 

1C34 

91 

THRU 

109 

1035 

1036 

1037 

1038 

1039 

1040 

1041 

ir.4 5 

lC4e 

104 7 

131 

THRU 

113 

1052 

1053 

1054 

1055 

1056 

1CS3 

1 122 

1 123 

1 124 

1125 

1 1 1 

THRU 

123 

1143 

1 144 

1145 

1 lUb 

2000 

2001 

2001 1 

20C12 

2002 

2003 

121 

THRU 

133 

2CP4 

20C5 

2006 

2007 

2006 

2C09 

2010 

201 1 

2C12 

201 3 

131 

THRU 

143 

20 14 

2015 

2016 

?U1 7 

?Clb 

2C19 

202C 

2021 

2u41 

?n« 3 
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141 

THRU 

150 

2044 

2045 

2046 

2047 

2048 

2049 

2050 2051 

2052 

2053 

151 

THRU 

16D 

2054 

2055 

2056 

2057 

2056 

2059 

2065 2066 

2067 

2070 

151 

THRU 

164 

2071 

2372 

2073 

2074 






CONDUCTOR ANALYSIS . . . 

, LINEAR = 

22, RtOIATlOM = 

1«»2, 

TOTAL = 

164, 

CONNECTIONS : 

716 



BCD JCONSUNTS D«TA 

NLOOP.SaOD ,nPLXCt, .OID.ARLXCA, .Oia.DAHPA ..B.DAHPO. .S 

END 


constants analysis... user = 

0, 

tOOEO =0 0 

0» TOTAL r 

0 



BCD 

3ARRAY DATA 







REM 

NODE ORDER FOR 

POyER 

DISTRIBUTION IS 

25, 

26, 

27, 

28, 

REM 

21, 22, 23, 24, 

101 , 

102, 103, 104, 

105, 106, 

107 

, 106 

, 109, 

REM 

no , 111, 112, 

113, 

114, 115, 116, 9 

* 10, 11, 

12, 

13, 

14, 15 

REM 

16, 29, 30, 31, 

32, 

17, 18, 19, 20 





REM 

HEAT INPUTS FOR 

A 1 

AMP BEAM CURRENT 

(BTU/HR ) 





1, 2.56, 2.56, 2.56, 2.56, 7.69, 7.69, 7.69, 7.69 
2.05, 2.C5, 2.05, 2.05, 6.20, 8.20, 8.20, 8.20 
10.25, 10.25, 10.25, 10.25, 102.5, 102.5, 102.5, 102.5 
19.73, 19.73, 19.73, 19.73, 2.05, 2.05, 2.05, 2.05 

16. m, 16.4 1, 16.41 , 16.41 , 2.56, 2.56, 2.56 , 2. 56, END 
REM HEAT INPUTS FOR A 2 AMP BEAM CURRENT (BTU/HR) 

2, 4.27, 4.27, 4.27, 4.27, 12.8, 12.8, 12.8, 12.8 

3.4, 3.4, 3.4, 3.4, 13.65, 13.65, 13.65, 13.65 

17.07, 17.07, 17.07, 17.07, 170.7, 170.7, 170.7, 170.7 
32.87, 32.87, 32.87, 32.87, 3.4, 3.4, 3.4, 3.4 

27.3, 27.3, 27.3, 27.3, 4.27, 4.27, 4.27, 4. 27, END 
REM HEAT INPUTS FOR A .5 AMP BEAM CURRENT (BTU/HR) 


4, 1.28, 

1.28, 

1.28, 

1.28, 

3.65, 

3.85, 

3.65 


3 . 

85 

1.03, 

1 .03, 

1.03, 

1.03, 

4.10, 

4.10, 

4.10 


4 . 

10 

5.12, 

5. 12, 

5.12, 

5.12, 

51.2, 

51.2, 

51.2 


51 

.2 

9.87, 

9.87, 

9.87, 

9.87, 

1 .U5, 

1.03, 

1.03 


1 . 

03 

8.20, 

8 . 20, 

3.20, 

8.20, 

1.23,. 

1.28, 

1.28 


1 . 

28 , END 

NODE ORDER FOR 

SOLAR 

FLUX 

IS 65, 

86, 87 

, 88 


25 

f 26, 


REM 29, 3C, 31, 32 

REM SOLAR INPUTS FOR 1 SUN INTENSITY 

REM THESE VALUES ARE DOUBLED FOR 2- SUNS INTENSITY 

3,30.0, 30. Q, 30.0, 30.0, 16.4, 18.4, 18.4, 18.4 

20. U, 20.0. 20.0, 20.0, END 

END 

ARRAY ANALYSIS... NUMBER OF ARRAYS = 4 TOTAL LENGTH : 136 


Execution Procedure 


BCD 3E XECUTION 

DIMENSION X(5000) F 

NDIM-50Ci0 F 

NTHrO F 

DIMENSION El(lOOC) ,E2llO0D),AREAl(lOD0),fiRF.A2(10O0),Fl?(10 0 0) F 

REM 

REM El ARRAY 
REM 

DATA (E1(I), 1 = 23, 164) / 5?«.6, 4U».15, .6, .775, .7, .675, 8*. 15, F 

1 . 1 75 , . 15 , . 15, .2 75 , 8 *. 1 5 , . 3 , 25m. 15 / E 

REM 

REM E2 ARRAY 
REM 
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DATA (: 2 m ,1=23,164) / 16*. 6Q, 76*. 15, F 

142*]., 3*1. ,5*1./ F 

5 EM 

PEM APEAl ARRAY (SQUARE INCHES) 

REM 

data ( area 1 ( I ) , I r23 , 1 64 ) / R.28, 12.2, 9.26, 12.2, 9. 28, 12.2, F 

1 2*9.26, 2*12.2, 7*9.23, 2*12.2, 2*9.28, 2*12.2, 2*9.28, 2*12.2, F 

2 2*9.28, 2*12.2, 2*9.28, 2*12.2, 2*9.28, 2*12.2, 2*9.28, 2*12.2, F 

3 2*9.28, 2*12.2, 2*9.26, 2*12.2, 9.28, F 

4 9.28, 2*12.2, 2*9.25, 2*12.2, 2*9.28, 12.2, 9.28, 12.2, F 

5 9*6.8, 3*12.7, 12.2, 19., 12.7, 2*13.8, 104.6, 2*12.1, 45.4, F 

6 2*3.98 , 2*4.4 , 12. 1 , 1 D4.6, 12.1, 104.6, '* F 

7 9.28, 2*12.2, °.29, 4*104.6, 4*12.1, 6.8,2*12.8, F 

8 48. ?, 4*33.1, 4*8.0, 28.8, F 

9 8*33.1, 8*8.0, 33.1, 104.6, F 

A 12.1, 4.4, 8.0, 29.8, 29.8, 2*4.4 / F 

REM 

REM AREA2 ARRAY (SOUARF INCHES) 

REM 

,i- \DAT'A ( AREA2 ri ) , I =23 , 1 64 ) / ,9.2e, 12.2 , 9.28 , 10*12.2 , 3*9.28 , F 

1 ' 12*6'. 77, 12*12 . 7, 12*34.7, 8*43.5 , 4*1 2.1 , F 

2 3*34.7, 3*12.7, 2*6.8, 3*34.7, 22.1, 33.3, 12.7, 2*12.1, 33.1, F 

3 8.0, 6.13, 252., 2*12.1, 2*8.0, 8.0, 33.1, 8.0, 33.1, 42*140., F 

9 8*140. / F 

REM 

REM internal VTEW FACTORS 
REM 

DATA (F 12( I ) , 1=23, 1 14 ) / .021, .027, .021, .023, .017, .016, .026, F 

1 .027, .026, .027, .0093, .021, .025, .0044, .013, .017, .0022, F 

2 .0035, .0059, .0088, .0084, .0105, .013, .0149, .0117, .0133, F 

3 .0148, .0156, .0296, .044, .071, .110, .022, .025, .026, .027, F 

4 ,0155, .0155, .0137, .010, ,437, .266, .153, .096, .043, .059, F 

5 .060, .054, .010, .029, .035, ,037, 1., 1., 1., .5, .0074, .0075, F 

6 .0073, .0038, . 50, .0098, , 00 3 8 , . 000 26 , .Q4S, .D3B, .048, F 

7 .166, .023, .006, .016, .016, .0785, .0425, .025, F 

a .108, .0041, 1., .22, .006, .302, .160, .23, 1., .55, .009, .39, F 

9 .003, .008, .020, .008, .020 / ^ 

REM / 

REM EXTERNAL VIEW FACTORS ' 

REM 

DATA (F 12( I ) ,1 = 1 15, 164 ) / .0891 , . 1095, . 137, . 174, E 

1 3*. 30, .30, 3*. 16, .16, .217, .5,1. ,.15, F 

2 3*.84, .84 , 3*. 71 , . 71 , 1 . , F 

3 .04 , .04, 3*. 04 , . 04 , 2*. 04 , . 031 , .031 , 3*. 031 , .031 , 2*. 031 , F 

. 4 .04 , .05 , .05 , .08 , . 24 , . 1 5, . 1 5 , . 39 , . 39/ F 

REM 

REM CALCULATE AND PRINT RADIATION CONDUCTOR VALUES ( B T U /H R . R ** 4 ) 

REM USING THE VALUES IN THE ABOVE ARRAYS 
REM 

DO 300 1=23,164 F 

IE(G(I).NE.O.) G.(I)=1. F 

300 CONTINUE F 

DO IDO 1=23,164 F 

G(I)=G(I)*1.19E-11*AREA1(I)*(1./((1./E1(I)-1.)+AREA1(I)/AREA2(I) F 

1*(1./E2(I)-1.)+1./F12(I))) F 

100 CONTINUE F 

WPITE(6,1D1) F 

101 FORMAT! 1 HI , 1 OX, 1 HI, 5X ,2HE1 , 5X,2HE 2, 5X, SHARE A1 ,5X, SHARE A2, ax, F 

13HF12,8X,1HG) F 
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WPrTEI6,lD2) (I,E](I),E?(I),«fit^l(I),AR£A?(I),n2(I) ,G(I),I:23,164 F 
1 ) F 

102 F0RKAT( 113 , 2F6. 2 ,F1 1 . 4 , riO.y ,2E14 .4 ) F 

REM CALCULATE STEADY STATE TEMPERATURES 
CINOSS 

PEM PRINT THRUSTER HEAT INPUTS 
QIPRNT 

END 

BCD 3VARIABLES 1 
REM 

REM STORE HEAT INPUTS FOR 2 AMP BEAM CURRENT 
' PEM USE SK ARAD ( A1+ 1 , ETC . ) FOR 1 AMP, 3 K A R A D ( A 4 ♦ 1 , E T C . ) FOR .S AMP 
REM 

BKARAD(A2+1, 025, 026, Q27f 02 fi, 021, 022,023,024 

QIC) 1,0 102,0103, 0104,0105,0 10 6,0107,01 0 8, 0109,01 in, Oil 1,01 12 
0113,0114,0115,0116,09,010,011,012,013,014,015,016,029,030 
C3 1 ,032,017,018,019,020) 

PEM 

REM ADO 3KARA0 ( A3+1,ETC) FOR SOLAR FLUX WHEN USED 

REM 

END 

BCD 3V ART ABLES 2 
END 

BCD 30UTPUT CALLS 

REM CONVERT temperatures FROM F TO K 


DO ion IR1,62 F 

Tm = 5./9.M ITU 1+43. 1-40. F 

T(I):T(I)+27i. F 

lOD CONTINUE F 

PEM PRINT TEMPERATURES IN K 
TPRINT 

REM CONVERT TEMPERATURES FROM K TO F 
00 200 IR1,62 F 

T(I):T(I)-273. F 

T ( I )Rl . 8M r T ( I ) + 411. ) -40. F 

2DD CONTINUE F 

REM PRINT TEMPERATURES IN F 
TPRINT 

END 


Radiation Properties and Calculated Radiation Conductors 


I 

El 

E2 

A R r A 1 

A R E A 2 

F 12 

G 

23 

.60 

.60 

9.28 0 0 

9. 2300 

.2100-01 

.2256-1 1 

24 

.60 

.6D 

12.2000 
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0.22791*02 

7 

111: 

0.227V1+C2 

T n?= 

5.22791*U2 

T 

1 1 J: 

7,5<a4?*?2 

T 

11“: 

7. 5264?. 3? 

? 

1 1 S: 

7.52662*02 

7 

110= 

’.02662*02 

7 

9r 

7.2676?*C2 

T JC = 

3,?b702*Ji 

I 

11 = 

3. 2675?. 07 

7 

1 ?: 

3.26757+C2 

7 

U: 

6,66061*3? 

7 

16 = 

•.66067*02 

7 

15 = 

6.«6«61»02 

T 16 = 

6,66061 *02 

I 

17 = 

5.00001*0? 

T 

1 *: 

“.00051*0? 

7 

19: 

0.00551*02 

7 

?n: 

0.oo00i*r? 

7 

21 = 

0.51714.R2 

7 22: 

5 , 5 1 ’ 1 6 ♦ 0 ? 

T 

2! = 

0.51?ia*O? 

T 

26: 

0.01216.3? 

7 

20 = 

3. 7030“*O2 

7 

?6 = 

3.70705*02 

7 

?7: 

3.?0’0*-*O: 

7 Pi: 

3.75355*P2 

T 

29: 

a.600Oe*O? 

7 

30: 

6.60008*02 

7 

31 = 

6.80008*0? 

7 

3?: 

•.60«?5*t2 

7 

33 = 

1 .6566' »n? 

7 36: 

1 .e0S60*L« 

I 

30: 

1.606NS+O? 

3 

56 = 

1 .60646,02 

1 

37: 

1.52882*02 

7 

5*: 

1 .026F?*C2 

7 

59 = 

1 .526t?*02 

7 4o: 

1.028b2*W 

T 

7a: 

7. 90067*02 

T 

75: 

?.85S67*.:2 

t 

76: 

2.90«0?*n2 

7 

77: 

2.e0P«7*O2 

7 

80 = 

1 » ji. a 6 a • Ci2 

T 06t 

l.r.''_’aeP*02 

I 

87: 

1. 30396*0? 

T 

S«: 

1.03868*32 

1 

66: 

2. 1 i3cr».o? 

7 

0Ci: 

-’.10003*32 

T 

*02 = 

-3.iOPtn*o:- 

7 *33: 

-3. lCCDO*a» 

1 

036: 

-3.10000*02 

7 

530: 

-!.13CCr*C2 












Lcopcr = 

90 fNGSU = 

. 1 7661 *00 












01 

1DI = 

J. 60330*03 

or 

I32r 

5.430'cn*00 

Cl 

103." 

3.6C0JfT*fl0 

01 

ir« = 

3.9Cn30*C0 

Of 

105 = 

1 . 360OO.C1 

n/ K6 = 

1 . 36500*01 

CI 

107: 

1 . 36000*01 

CI 

108S 

1 . 360or*oi 

0! 

ICV = 

1 . 70700*01 

01 

110= 

1.707C0.CI 

OI 

111 = 

1 .7070’?. 01 

01 11? = 

1.73700*01 

CI 

113 = 

1.7ti70C*0? 

01 

1165 

1 . 70700*02 

CI 

1 10 = 

1.7C70C**J? 

CI 

116 = 

1 .70700*02 

OI 

9: 

?.?87U*>*0I 

CI 10 = 

;.?870P*01 

CI 

11 = 

3. 26700*01 

01 

125 

3.?870n*Sl 

01 

13 = 

3.«Ciaoc*oo 

01 

19 = 

3.9cncc»*tjo 

OI 

15 = 

3.6000'C*Ci. 

OI 16 = 

3.atnnco*0u 

CI 

17 = 

4.27000*03 

01 

18:: 

8.273cn*00 

01 

19: 

6.2700C*CU 

Cl 

20= 

9, 27000*00 

OI 

21 : 

1 .?onc~*ci 

CI ?? = 

l.26000*Cl 

01 

73 = 

1.2fe00C«Cl 

01 

26; 

I .28000*01 

CI 

20: 

6,27000*00 

31 

26 = 

9.?7P00*C6 

CI 

27 = 

6.27PC0*Cj 

CI 28 = 

4.270or*CO 

01 

29 = 

2. 733:r*01 

01 

30? 

?,73DCn*01 

;i 

31 = 

2. 7300n*01 

Cl 

*?: 

?.73oro*oi 

01 

53 = 

O.RCDCr 

01 36: 

z . ci-ncn 

CI 

*5 = 

C. COQOO 

01 

’6? 

0. 00000 

CI 

37 = 

O.OCPOD 

01 

’8 = 

o.popcn 

OI 

39 = 

n.Punar. 

01 6C = 

c.ooncc 

01 

76 = 

o.ocoon 

01 

70? 

c.oocon 

01 

76 = 

o.oucor 

Oi 

77= 

''.OOPCO 

01 

8*: 

n.ocoo'’ 

01 5t = 

L.POPCG 

01 

67 = 

o.cuooo 

CI 

68? 

n.oooun 
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TABLE I. - DESCRIPTION OF THE ANALYTIC NODAL NETWORK 


Node 

Location 

Material 

101 - 104 

Anode, rear 

6A1-4V titanium 

105 - 108 

Anode, near rear 



109 - 112 

Anode, near front 



113 - 116 

Anode, front 



9 - 12 

Engine body, rear and structure 



13 - 16 

Engine body, front and ring 

6A1-4V titanium and 1010 carbon steel 

17 - 20 

Pole 

1010 carbon steel 

21 - 24 

Screen grid 

Molybdenum 

25 - 28 

Accelerator grid 

Molybdenum 

29 - 32 

Baseplate 

6A1-4V titanium 

33 - 36 

Ground screen, rear 

6061-T6 aluminum 

37 - 40 

Ground, screen, front 

6061-T6 aluminum 

46 

Rear shield 

6061-T6 aluminum 

74 - 77 

Grid mounting assembly 

304 stainless steel 

85 - 88 

Mask 

Pure titanium 

501 - 504 

Chamber or shield boundary 



505 

Chamber boundary 




TABLE II. - ASSUMED PHYSICAL PROPERTIES OF 


ION THRUSTER MATERIALS 


Material 

Density, 

g/cm^ 

Specific heat 
at 573 K, 
J/(g)(K) 

Thermal con- 
ductivity 
at 573 K, 
W/(cm)CK) 

6061-T6 aluminum 

2.77 

0.837 

1.80 

Pure titanium 

4.43 

.628 

.20 

Carbon steel 

7.81 

.544 

.60 

304 stainless steel 

7.92 

.523 

.20 

Molybdenum 

10.19 

.837 

1.20 

Tantalum 

16.16 

. 146 

.60 

Mercury 

13.56 

. 126 

. 10 

Tungsten 

19.38 

. 146 

1.50 

Alumina (Al20g) 
(WESGO Al-300) 

3.79 

.837 

.17 

Kovar 

8.36 

. 439 

.15 

6A1-4V titanium 

4.43 

.628 

.10 





TABLE m. - ASSUMED SELF-HEATING 


POWER DISTRIBUTIONS 


1 

Component 

Ion beam current, A 

Percent of 
total power 

0.5 

1.0 

1 

2.0 

Self- heating power, W 

Accelerator grid 

1.5 

3.0 

5 

1.47 

Screen grid 

4.5 

9.0 

15 

4.42 

Anode, rear section 

1.2 

2.4 

4 

1.18 

Anode, middle sections 

10.8 

21.6 

36. 0 

10.60 

Anode, front section 

60.0 

120.0 

200.0 

58.91 

Engine body, rear 

11.6 

23.7 

38.5 

11.34 

Elaine body, front 

1.2 

2.4 

4. 0 

1. 18 

Baseplate 

9.6 

19.2 

32. 0 

9.43 

Pole 

1.5 

3.0 

5. 0 

1.47 

Total 

101.9 

203.8 

339. 5 

100.00 

Discharge power 

111 

222 

370 





TABLE IV. - COMPARISON OF EXPERIMENTAL AND ANALYTIC TEMPERATURES^ FOR THE 


CALIBRATION TESTS WITH NO SHIELD OR SOLAR INCIDENCE 


Component 

Ion beam current, A 


2.0 

1.0 

0.5 


T 

e' 

T 

a* 

Error, 

T 

e’ 

T 


Error, 

T 

T 


Error, 


K 

K 

- Ta) 

K. 

K 

(Tp - T„) 

K 

K 

(T. - T„) 





° ^ 100, 




° ^ 100, 




— — 100, 









T 




T„ 





e 




e 




e 





percent 




percent 




percent 

Baseplate 

513 

514 

-0. 19 

467 

455 

2.57 

408 

384 

5.88 

Engine body, rear - 













Under neutralizer housing 

... 

438 



388 


... 

330 

— 

.. 90° from neutraUzer 

439 



.23 

402 



3.48 

348 



5. 17 

180° from neutralizer 

448 



2.23 

408 



4.90 

354 



6.78 

270° from neutralizer 

439 



.23 

402 



'3.48 

349 



5.44 

Engine body, front - ' 













Under neutralizer housing^ 

516 

504 

2.33 

465 

446 

'4.09 

396 

376 

5.05 

90° from neutralizer 

493 



-2.23 

446 



0 

378 



.53’ 

180° from neutralizer^ 

511 



1.37 

458 



2.62 

387 



2.84 

270° from neutralizer 

501 



• -.60‘ 

450 



.89 

378 



.53 

Anode, rear 

517 

523 

-1.16 

466 

462 

.86 

401 

389 

2.99 

Anode, front 

674 

673 

.15 

610 

•591 

3.11 

518 

494 

4.63 

Rear shield^ 

373 

373 

0 

356 

356 

0 

317 

317 

o' 

Screen mask, 90° from neutralizer 

332 

312 

6.02 

310 

281 

9.35 

263 

243 

7. 60 , 

Ground screen, center 













90° from neutralizer 

320 

349 

-9.06. 

303 

307 

-.1.32 

267 

259 

3. 00 

180° from neutralizer 

332 

349 

-5.12 

311 

307 

,1.29 

273 

259 

5. 13 

270° from neutralizer 

329 

349 

-6.08 

311 



307 

1.29 

272 

259 

4. 08 

Average absolute error 


2. 64 


2.80 


4. 26 


^Tg = experimental temperature; T^ = analytic temperaWe. 
^Thermocouples not on engine body, but on ring near grids. 
*^Set to average experimental values. 


TABLE V. - SOLAR HEATING DISTRIBUTION 


[Q^^js = NSAofg where N is the number of suns incident, S is the 

solar constant (= 0. 135 W/cm^), A is the absorbing area,' cm^, 

and or„ is the effective solar absorptivity.] 
s 


Component 

Solar Flux 

“s 

1 

sun 

2 suns 

Incident, 

W 

Absorbed, 

W 

Incident, 

W 

Absorbed, 

W 

Mask 

76.0 

35.2 

152.0 

70.4 

0.46 

Accelerator grid 

44.9 

21.6 

89.8 

43.2 

.48 

Baseplate 

44.9 

23.4 

89.8 

46.8 

.52 
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TABLE VI. - COMPARISON OF EXPERIMENTAL AND ANALYTIC TEMPERATURES^ FOR THE 
CALIBRATION TESTS WITH NO SHIELD AND 2 SUNS SOLAR INCIDENCE 


Component 

Ion beam current, A 

2.0 

1.0 

0.5 

Te> 

K 

Ta* 

K 

Error, 

(T^ - T„) 

® ^ 100 

Te’ 

K 

T 

a’ 

K 

Error, 

(T„ - T ) 

® ** 1^00 

T 

e’ 

K 

Ta> 

K 

Error, 

6 ® 100 

Te 

percent 

Te 

percent 

' percent 

Baseplate 

549 

553 

-0.73 

508 

507 

0.20 

469 

461 

1.74 

Engine body, rear - 













Under neutraUzer housing 


455 



412 


— 

368 


90° from neutralizer 

466 



2.36 

431 



4.41 

394 



6.60 

180° from neutralizer 

477 



4.61 

439 



6. 15 

402 



8.46 

270° from neutralizer 

467 



2.57 

431 



4.41 

394 



6.60 

Engine body, front - 













Under neutralizer housing'* 

548 

529 

3.47 

503 

480 

4.57 

456 

430 

5.70 

90° from neutralizer 

527 



-.38 

482 



.41 

434 



.92 

180° from neutralizer** 

540 



2.04 

494 



2.83 

446 



3.59 

270° from neutralizer 

533 



.75 

486 



1.23 

436 



1.38 

Anode, rear 

550 

550 

0 

507 

498 

1.78 

463 

445 

3.89 

Anode, front 

698 

686 

1.72 

636 

609 

4.25 

554 

524 

5.42 

Rear shield** 

394 

394 

0 

371 

371 

0 

347 

347 

0 

Screen mask, 90° from neutralizer 

468 

467 

.21 

461 

458 

.65 

450 

449 

.22 

Ground screen center 













90° from neutralizer 

341 

394 

-15.54 

318 

369 

-16.04 

297 

347 

-16.84 

180° from neutralizer 

352 

394 

-11.93 

330 

369 

-11.82 

307 

347 

-13.03 

270° from neutralizer 

352 

394 

-11.93 

330 

369 

-11.82 

308 

347 

-12.66 

Average absolute error 


4.16 


5.04 


6.22 


^T = experimental temperature; T = analytic temperature. 

C a. 

Thermocouples not on engine body, but on ring near grids. 
*^Set to average experimental values. 
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TABLE Vn. - COMPARISON OF EXPERIMENTAL AND ANALYTIC TEMPERATURES'^ FOR MULTIPLE THRUSTER SIMULATION TESTS 



CO CO CD CO 
CO CO CO 

m I/O m 


^ 00 
C>J ^ CO 
LA t /0 t /0 


t-COOCO^OO^CvJ 

rrcDC-CDC-OiCDO^ 

t/Ol/OL/OlOlfOCD^^ 


i/ 0 *ACOC'JC-*H^^ 
tDC— t-C^C—i^CDO 
l/Jl/OlAl/Ol/OC“^l/0 


^ CO 
05 O cs 
^ m ^ 


UO CO CO 
o O 05 
I/O I/O 


O 04 I/O CM 
C- CM 04 C4 
V I/O I/O I/O 


CO CO C4 
*-( CM CM 
I/O I/O I/O 


lAl/OC“l/0 OCM^O 
eoi/ouoi/ocoo>L/o^ 
L/0l/0«/0l/0l/0CO'0*^ 


fHOcot-t-^-ii-ieo 

l/5CDCDlfOUO«-^UOl/0 

i/ 5 l/Ol/Ol/ 5 l/OC-^^ 


CO CM 00 
C- CJ5 C- 
^ ^ ^ 


05 *H 05 
CO 05 
^ ^ ^ 


05 05 CM 05 

lA O »-H O 
^ m I/O t /0 


I 05 O CD 

C 05 «-t O 

C ^ t/0 I/O 


«>iOeMOCD^OcD 

l/Oi/Ol/Oi/Ol/OCD^CO 


^eot*»-' 0 oc 000 

CO^^^^ 05 ^C 0 

i/oi/oi/oi/oi/oco^eo 


CD CM O O CO CO 

CO O CM CM CO *H 


t- 

I/O t- I/O 
^ 


t“ 05 t“ 
t- C- CO 
^ ^ ^ 


CO t* 05 ^ 

m ^ ^ I/O 

I/O L/o tro t/o 


CM CO CO 05 
CO O f 05 
LfO t- ■M* 


CMCDt/Ol/OCOi^COCO 
OOOOCOCOCO^C— 05 
l/OC/Ot/OI/Ol/OC~^^ 


CM CM CM 
CM CM CM 
I/O LfO LfO 


O O CM 
CM 1-H ^ 
I/O I/O I/O 



w -M* O 
CM W 1-* CM 
C/0 L/O I/O I/O 


CO CM CO O C- 

'M* — ^ I/O CO m c> 

I/O I/O CO ^ CO 


CO’— 'CMCOO^OC- 

i/ocoi/oiroi/ocMi/oco 

l/OI/Ol/Ol/OI/OC-^CO 


CO CO CO 
CO 00 CO 
^ ^ ^ 


uro CO '-H 
CO c- op 
^ 


flj *3 

§ I 

fcl P 

•O g s 

I 5 : 2 

i s |o ' 

0 ) .5 c ^ 

ffl gi D a. 

O W 


a) 5 
s « 2 


•sis 


® -S ' 

^ P 05 

u 


V 05 
05 N tM 

.2 sa ;3 

^ ci oj 

S B B 

XI 9 9 


t. "S “ 

o 3 

2 i ® 


® 5 " 
2 ® ” 
S S I 


S o o o 
3 o CO c- 

p 05 CM 
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T = experimental temperatures; T = analytical temperatures. 

I ^ ** 

Thermocouples not on engine body, but on ring near grids. 

'Set to average experimental values. 












































TABLE VIII. - COMPARISON OF EXPERIMENTAL AND ANALYTIC TEMPERATURES^ FOR MULTIPLE THRUSTER SIMULATION TESTS 



^ t- ^ 
CO CO CO CO 
^ ^ TT 


t- Oi <0 

00 o o 

Tl* TT 


coc“ 0 >t-cge>J^O 

ost-HWi-^eqc'jcoco 


caOiC'jirtCMO'^co 

egc'3oocMco4r?coco 


CD ^ 

mom 

^ TJ4 


m m m 
CD CO m 

T}< ^ ^ 


1 -^ cvi cr- o tr- 

O CnJ CD c- CD 

m Tl< Tj« tT TT 


<Ti I m m 1-H 

O I CD C- t- 

m I ^ ^ ^ 


i-( CO O CO 04 CD vH ■ 

coo>ooiOiHi-Hesj 
.^^m^mcD^^ 


t-COCOOt-COCDiH 

a^OOooco<- 4 co 

•crmmmmcD^^ 


1-4 04 

CO CO 
TP 


04 ^ tP 

tP tP CO 
tP ^ 


m 04 ^ CO 
O ^ TP ^ 
^ TP ^ 


m ^ »H 
^ m m 

TP TP ^ 


o^comcocooooo 

mc-c«i>c-ocoo 4 

Tp^Tp^TpCDCOCO 


CDOCDC 0 ^ 0 - 0 ) 0 ) 

c-cocot-cotMcoeo 

^^^^^CDCOCO 


O w o 
O) o o> 
eo ^ CO 


CO O 05 c- 
^ th O w 

m m m m 


T-i 04 m O 


04 CO 05 ^ 05 

CO ' » CO 04 m CO 

m m CD 


050 »MO>^ 05 ^CD 

^m^co'^cDmco 

mmmmmcD^^ 


r- «-i 00 CD 
r- p- CD c« 

Tp TP TP TP 


m 00 »-< CD 

05 ^ o 04 m 

^ CD CO 


mcDcoir-coT-ii-<05 

00050505 ^ 04 m 

mm^^^co^co 


o- r- c- 

Tp TP tp 
^ TP TP 


CO 

' ■ i 

X 05 

. ^ 

a> t — 4 

^ c 4 ^ 

t S b 

L rt 3 
0 > 

■a I S 
Si 
tt) .S ? ^ 

« g) P 05 


J 3 

(h U 
(h V Q> . 

05 N N 

rt ctf 

^ b b 

^33 
3 0> 0) 

0 ) C 3 


S ' M U 

^ ^ 05 05 

O S ® « N 

^ -a 2 P p 

o_ o g is is 

S is 3 3 


LBEijc'io’agzga 

6 ooSo 2 “S 3 Soo 

o ” S §o 'fe'o 

OOOt-^^ ^ 05 Socot- 

T < K w O 
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^T- = experimental temperature: = analytic temperature. 

»i a 

Thermocouples not on engine body, but on ring near grids. 
Set to average experimental values. 











































TABLE DC. - COMPARISON OF EXPERIMENTAL AND ANALYTIC TEMPERATURES® FOR MULTIPLE THRUSTER SIMULATION TESTS 
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Tg = experimental temperature; = analytic temperature. 
’Thermocouples not on engine body, but on ring near grids. 
Set to average experimental values. 












Figure 3. - Schematic of 30-centimeter 
diameter ion thruster with 360° 
azimuthal shield. 


(-Ground screen 



Figure 4. - Schematic of 30-centimeter 
diameter ion thruster with 27(P 
azimuthal shield. 


312 

(332>-\ 



Figure 5. - Analytic and experimental (in parentheses) temperatures (kelvins) for 2-ampere 
beam current with no shield or solar incidence. 



(387)^ 
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Figure 8. - Analytic and experimental (in parentheses) temperatures (kelvins) for 2-ampere Figure 9. - Analytic and experimental (in parentheses) temperatures (Kelvins) for 2-ampere 

beam current with 360° azimuthal shield at 525 K and 2 suns solar incidence. beam current with 27CP azimuthal shield at 483 K and no solar incidence. 




Figure 10. - Analytic and experimental (in parentheses) temperatures (Kelvins) for 2-ampere 
beam current with 270° azimuthal shield at 498 K and 1 sun solar incidence. 
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Figure 11. - Analytic and experimental (in parentheses) temperatures (Kelvins) for 2-ampere 
beam current with 270P azimuthal shield at 464 K and 2 suns solar incidence. 
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to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 




